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We present theoretical calculations for polarization and ellipticity of high-order harmonics from 
aligned N2, CO2, and O2 molecules generated by linearly polarized lasers. Within the rescattering 
model, the two polarization amplitudes of the harmonics are determined by the photo-recombination 
amplitudes for photons emitted parallel and perpendicular to the direction of the same returning 
electron wave packet. Our results show clear species-dependent polarization states, in excellent 
agreement with experiments. We further note that the measured polarization ellipse of the harmonic 
furnishes the needed parameters for a "complete" experiment in molecules. 
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High-order harmonic generation (HHG) is one of the 
most important nonlinear processes that occur when 
atoms or molecules are placed in an intense laser field [1[ . 
Today these high harmonics are used as convenient lab- 
oratory XUV or soft X-ray light sources, as well as the 
sources of single attosecond pulses or attosecond pulse 
trains 0-0] • High harmonics are emitted when laser- 
induced continuum electrons recombine with the target 
ions. Since photo-recombination is a time-reversed pro- 
cess of photoionization (PI), study of HHG from molecu- 
lar targets offers alternative means for probing molecular 
structure that have been traditionally carried out using 
PI at synchrotron radiation facilities. Gaseous molecules 
can be given a periodic transient alignment by a weak 
short laser pulse [5]. By studying HHG generated from 
such aligned molecules, information such as molecular 
frame photoelectron angular distributions (MFPAD) for 
PI from valence orbitals of molecules can be inferred. The 
goal of a "complete experiment" is to determine ampli- 
tudes and phases of all dipole matrix elements. For lin- 
ear molecules, this may be achieved if measurements of 
MFPAD are carried out using elliptically polarized lights 
For photo-recombination, this means that one may 
obtain equivalent information by examining the elliptical 
polarization of HHG from aligned molecules. 

Clearly if the gas is isotropically distributed, as for 
atomic or unaligned molecular targets, due to the sym- 
metry the emitted harmonics are polarized parallel to 
the polarization of the driving linearly polarized laser. 
For aligned molecules, a harmonic component perpen- 
dicular to the laser polarization is expected to be present 
in general [8j. This requires that experiments be car- 
ried out with a good level of molecular alignment in or- 
der to observe a significant amount of the perpendicular 
harmonic component. It is therefore not surprising that 
polarization measurements for emitted harmonics were 
reported only very recently (9l4l3| . All these experiments 
were carried out within the pump-probe scheme, where 
a relatively weak, short laser pulse is used to impulsively 



align molecules along its polarization direction, and after 
some delay time, a second laser pulse is used to generate 
high-order harmonics. We note that the commonly used 
strong- field approximation (SFA) cannot be employed to 
interpret such experiments since it predicts little or no 
ellipticity for emitted harmonics 0, [l4| . 

In this Letter we report theoretical results for polar- 
ization and ellipticity of HHG from aligned N2, O2, and 
CO2 molecules. Our results show very good agreement 
with experimental measurements 0, [To|, El for harmonic 
orientation angles and the reported large ellipticity for N2 
[13| . The calculations were carried out using the quanti- 
tative rescattering theory (QRS) [l5l-[l9| where the com- 
plex induced dipole responsible for harmonic emission is 
represented as a product of a returning electron wave 
packet and the laser- free photo-recombination transition 
dipole, 



D h± (u,0)=W(E k ,0)d\ U± (w,0). 



(1) 



Here $ is the angle between the molecular axis and the 
(probe) laser polarization direction, Ek is the "incident" 
energy of the returning electron, and uj = I p + Ek is the 
emitted photon energy, with I p being the ionization po- 
tential. The returning electron can recombine with the 
parent ion to emit a photon with polarization in the par- 
allel or perpendicular direction to its motion, resulting in 
the two polarization components of the emitted harmon- 
ics. Both of these complex transition dipoles d^± are 
obtained from state-of-the-art molecular photoionization 
code for each fixed- in- space molecule. Note that 

we only need to consider the harmonic components on the 
plane perpendicular to the propagation direction of the 
driving laser, since only the harmonic emission propagat- 
ing along this direction can be efficiently phase matched. 
As for the returning electron wave packet, we extract it 
from the SFA [lfj. Eq. (1) thus shows that the ampli- 
tude and phase of the transition dipoles can be probed 
by studying HHG. 

To compare with experiments, induced dipoles 
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FIG. 1: (Color online) Fixed- in-space molecular photoioniza- 
tion differential cross sections, corresponding to the parallel 
and perpendicular components of emitted harmonics H17 and 
H23 for N 2 (a), C0 2 (b), and 2 (c). 



from the fixed- in- space molecules are coher- 
ently convoluted with the molecular alignment distribu- 
tions [sl EH, [Hj. We note that this alignment "phase- 
matching" tends to favor the parallel component. In our 
simulations, the alignment distribution is obtained from 
numerical solution of the time-dependent Schrodinger 
equation within the linear rotor model for each molec- 
ular species 0, [l6[. We use a 120 fs pump laser pulse 
with an intensity of 3 x 10 13 W/cm 2 , and a 30 fs probe 
laser pulse with an intensity of 2 x 10 14 W/cm 2 . Both 
pump and probe lasers are of 800 nm wavelength. The 
rotational temperature is assumed to be 100 K. These pa- 
rameters were chosen to closely match the experimental 
conditions of Zhou et al [13[ . We vary the angle between 
pump and probe polarizations and use the half-revival 
time delay for N2 and O2, and 3/4-revival for CO2. 

To understand the experimental measurements, we 
show in Fig. 1 fixed- in- space PI differential cross sections 
for the three species, corresponding to the parallel and 
perpendicular components of emitted harmonics HI 7 and 
H23. Note that these cross sections are proportional to 
|d||| 2 and |<i±| 2 , respectively. In general, the perpendicu- 
lar components are smaller than the parallel components 
for all three targets. Due to symmetry, the perpendicular 
harmonic component will vanish after averaging over the 
alignment distribution if the angle between pump and 
probe polarizations is 0° or 90° (see Fig. 2 below). 
From Fig. 1 we note that at intermediate angles the cross 
section for the perpendicular component is quite compa- 
rable to the parallel one for N2 and CO2. For O2 the 
perpendicular is always much smaller than the parallel 
one. Therefore one can expect a small intensity ratio be- 
tween perpendicular and parallel components for O2, but 



FIG. 2: (Color online) Experimental (left panels) and theoret- 
ical (right panels) intensity ratios between perpendicular and 
parallel component of harmonic fields, as a function of align- 
ment angle between pump and probe polarization directions, 
for N 2 (top panels) and CO2 (bottom panels). Experimental 
results are taken from Zhou et al 11311 . 



a larger ratio for both N2 and CO2. 

Next we show in Fig. 2(b) and 2(d) theoretical inten- 
sity ratio j^- between perpendicular and parallel compo- 
nents for harmonic orders from H17 to H23, as a function 
of alignment angle between pump and probe polariza- 
tions. Experimental results by Zhou et al [l3[ are shown 
in Fig. 2(a) and 2(c) for comparison. For N2, the theoret- 
ical intensity ratio reaches 0.16 at the peak near 55° for 
all the harmonics from H17 to H23 [Fig. 2(b)]. The mea- 
surement [13| shows a very similar shape with peaks near 
55° as well, but the magnitude increasing with harmonic 
orders from 0.1 to 0.22. For C0 2 [Fig. 2(d)], the shape of 
the intensity ratio from theory changes slightly, with the 
peak now at about 40°. This is also in good agreement 
with experiment. As for the magnitude, the theoretical 
intensity ratio is about a factor of three larger than in the 
experiment. This discrepancy could be partly due to the 
fact that near the minimum of the parallel component 
[see Fig. 1(b)] the theoretical transition dipole is not ac- 
curate enough or that inner molecular orbitals may con- 
tribute [23|,[24|. For O2, we found that the intensity ratio 
is very small, as expected, with the biggest intensity ra- 
tio of about 1% near 35°. We comment that for all three 
targets, the ratio goes to zero if pump and probe polar- 
izations are parallel or perpendicular. As stated earlier, 
this is expected from symmetry consideration. 

Let S be the phase difference between perpendicu- 
lar and parallel components of the harmonic field and 
tan(7) = \J~^- As S 7^ or tt in general, the emitted 
harmonic is elliptically polarized. To characterize the po- 
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FIG. 3: (Color online) Experimental (left panels) and the- 
oretical (right panels) orientation angle (j) (in degrees) as a 
function of alignment angle between pump and probe polar- 
ization directions and harmonic order for N2 (top panels), 
CO2 (middle panels). Experimental results are taken from 
Zhou et al 13]. Theoretical result for O2 is also showed (e). 

larization ellipse, we define [25j the orientation angle (f) 
of the ellipse and the ellipticity e = tan(x) by 

tan(20) = tan(27) cos(£), (2) 
sin(2x) = sin(27) sin(5). (3) 

Our results for the orientation angle <fi are shown in 
Fig. 3(b) and 3(d), as a function of alignment angle be- 
tween pump and probe polarizations and harmonic or- 
der. Experimental results by Zhou et al [13] are also 
shown (left panels) for comparison. The theoretical data 
are anti-symmetric with respect to the sign change in the 
pump-probe angle 0, so in the following we just focus on 
the positive 6. The experimental data are less symmet- 
ric. Here the positive (or negative) angle corresponds to 
the case when the major axis of the ellipse rotates in the 
same (or opposite) direction as the molecular axis. 

The most noticeable feature for N2 is the sign change 
in the orientation angle as a function of harmonic order 
near HI 9. The orientation angle of about 20° at H13, de- 
creases smoothly with harmonic order, and reaches —20° 
at H29. This is in excellent agreements with the measure- 
ments by Zhou et al, shown in Fig. 3(a), as well as with 
Levesque et al [9]. Zhou et al [13] found that the sign of 
the orientation angle changes near H19, while Levesque 
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FIG. 4: (Color online) Experimental (left panels) and theoret- 
ical (right panels) ellipticity (top panels) and phase difference 
(bottom panels) between the two polarization components 
from N2 for angles between pump and probe polarizations 
— 40°, 50°, 60°. Experimental results are taken from Zhou 
et al . 

et al [9] found the change near H21, independent of the 
pump-probe polarization angle. We comment that cal- 
culations based on the SFA do not lead to a satisfactory 
agreement with experiments [9]. For CO 2, the theoretical 
orientation angles are negative (for positive 0) for all the 
considered harmonics. This is in good agreements with 
Zhou et al, shown in Fig. 3(c), and Levesque et al Q. 
For O2, on the other hand, the orientation angle remains 
positive for the harmonic range shown in the figure. Its 
magnitude is also much smaller, reaching about 6° near 
H19-H21 for ~ 30° - 40°. This behavior agrees well 
with Levesque et al [9j]. Again, these small orientation 
angles are due to the small intensity ratios, which in turn 
is related to the smallness of the PI cross sections for the 
perpendicular component, as compared to the parallel 
component, for intermediate angles in O2 [see Fig. 1(c)]. 

We also compare ellipticity e and phase difference S vs 
harmonics order at fixed pump and probe angles = 40°, 
50°, and 60°. Fig. 4 shows that the theoretical results 
for N2 are in good agreement with the experimental data 
of Zhou et al [13]. In particular, the theory predicts a 
large ellipticity up to e ~ 0.4 near H21, in agreement 
with experiment. We comment that a recent calculation 
based on an extended stationary-phase SFA by Etches 
et al [14] showed very weak ellipticity of about 0.02 only. 
The phase difference, shown in Fig. 4(d), increases nearly 
linearly with harmonic order, from 0.17T at H13 to 0.87T at 
H27, but is nearly independent of alignment angle. Note 
that the phase difference is nearly it/ 2 at HI 9. This is 
exactly the harmonic order, where the orientation angle 
changes its sign, see Eq. (2) and Fig. 3(b). This behav- 
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ior is in excellent agreement with experiment shown in 
Fig. 4(c). For C0 2 the ellipticity from the QRS is slightly 
smaller than that of N2, while the measurements by Zhou 
et al showed a value of less than 0.1. This discrepancy 
is consistent with the larger errors we found for the CO2 
intensity ratio, but the reason remains largely unclear at 
this moment. We further note that the calculation for 
CO2 by Smirnova el al [23| showed an ellipticity of 0.1 at 
H29, which increases up to about 0.4 at harmonics H37- 
H43. In their simulation, contributions from two lower 
molecular orbitals were also included. For completeness 
we note that the QRS predicts an ellipticity for O2 of less 
than 5% under the same experimental conditions. 

In general, experimental HHG spectra include the ef- 
fect of macroscopic propagation in the medium (26| . 
However, under typical experimental conditions, we can 
show that macroscopic propagation will affect both har- 
monic components in the same way. Indeed, the propa- 
gation equation for each harmonic component E a (with 
a =|| or _L ) ca n be written under the paraxial approxi- 
mation as [26|, 127 ] 



2iuj dE a (r, z, lj, 0) 
c dz 



oc (D a (r,z,u))o,{4) 



where (-D|| } j_(r, z,u))q is the nonlinear polarization, aver- 
aged over the molecular alignment distribution for a fixed 
pump-probe angle 0. Here we assume that absorption 
and free-electron dispersion are negligible. In a typical 
gas jet experiment, the aligning laser is much less intense 
and more loosely focused than the probe laser. There- 
fore we can assume that the aligning laser is uniform 
in the gas jet, which is typically of about 1 mm thick. 
We found that for a fixed {oj,0} the intensity ratio and 
phase difference between the two components (D\\(u))q 
and (D±(uj))o change less than 10% as probe laser in- 
tensity changes from 1.5 x 10 14 to 2.5 x 10 14 W/cm 2 . 
In other words, the ratio R = \(D±) / (D\\)\ and phase 
difference are nearly independent of the spatial coordi- 
nates {r, z} in the gas jet. From Eq. (4), it follows that 
the ratio \E±/E\\\ = R. The same arguments also hold 
for the phase difference between E± and En. This im- 
plies that the results presented in this paper should be 
nearly unchanged even if the macroscopic propagation is 
carried out. Our results are still dependent on the de- 
gree of molecular alignment, which is controlled by the 
pump pulse. Therefore polarization resolved HHG mea- 
surements allow us to directly extract single- molecule fea- 
tures (up to averaging over the alignment distribution) 
without much influence of the details of the macroscopic 
phase-matching conditions. 

In conclusion, we have shown that the quantitative 
rescattering theory can be extended to calculate polariza- 
tion and ellipticity of high-order harmonics from aligned 
molecules in intense laser fields. Theoretical results are 
compared to experimental measurements side by side and 
good agreement has been found. The interaction of light 



with molecules is governed by the dipole transition ma- 
trix elements. This dipole interaction has been tradi- 
tionally probed using photoionization, but can similarly 
be probed by studying HHG. While photoionization has 
the advantage of achieving higher energy resolution to 
reveal many-electron dynamics, HHG has the advantage 
of surveying a broader photon energy range coherently 
in one single experiment, thus revealing the global prop- 
erty of the molecule. Since the phases of the harmonics 
can be conveniently measured experimentally, HHG also 
has the advantage of revealing directly the phases of the 
transition dipoles. 
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